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ABSTRACT 



We interpret the rapid correlated UV/optical/ X-ray variability of 
XTE Jl 118+480 as a signature of the coupling between the X-ray corona and a 
jet emitting synchrotron radiation in the optical band. We propose a scenario in 
which the jet and the X-ray corona are fed by the same energy reservoir where 
large amounts of accretion power are stored before being channelled into cither 
the jet or the high energy radiation. This time dependent model reproduces the 
main features of the rapid multi-wavelength variability of XTE J1118+480. As- 
suming that the energy is stored in the form of magnetic field, we find that 
the required values of the model parameters are compatible with both a patchy 
corona atop a cold accretion disc and a hot thick inner disc geometry. The range 
of variability timescales for the X-ray emitting plasma are consistent with the 
dynamical times of an accretion flow between 10 and 100 Schwarzschild radii. 
On the other hand, the derived range of timescales associated with the dissipa- 
tion in the jet extends to timescales more than 10 times larger, confirming the 
suggestion that the generation of a powerful outflow requires large scale coherent 
poloidal field structures. A strong requirement of the model is that the total jet 
power should be at least a few times larger than the observed X-ray luminosity, 
implying a radiative efficiency for the jet e- 3 < 3 x 10~ 3 . This would be consistent 
with the overall low radiative efficiency of the source. We present independent 
arguments showing that the jet probably dominates the energetic output of all 
accreting black holes in the low-hard state. 

Key words: accretion, accretion discs - black hole physics - magnetic field - 
radiation mechanisms: non-thermal - star: XTE Jl 118+480 - X-rays: binaries 



1 INTRODUCTION 

The high energy spectrum (> 1 keV) of accreting stellar 
mass black holes in the low/hard state can be roughly 
described by a power-law with photon index r ~ 1.4 — 2, 
and a nearly exponential cut-off at a characteristic energy 
E c of a few hundred keV (see e.g. Tanaka & Lewin 1995; 
Gierlinski et al. 1997; McClintock & Remillard 2004). 
Such a spectrum is generally interpreted as due to ther- 
mal Comptonisation in a plasma with electron tempera- 
ture kT e ~ 100 keV and Thomson optical depth t ~ 1 
(see e.g. Poutanen 1998). There are two possible explana- 
tions for the presence of this very hot plasma (often called 
corona). It could be either a geometrically thick, optically 
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thin innermost part of the accretion flow (Shapiro, Light- 
man & Eardley 1976; Narayan & Yi 1994) or a collection 
of small scale active regions located atop a cold, geomet- 
rically thin and optically thick accretion disc (Haardt, 
Maraschi & Ghisellini 1994), possibly powered by mag- 
netic reconnection. 

When the X-ray luminosity increases above a few 
percent of the Eddington luminosity (I/Edd), accreting 
black holes are observed to switch from the low/hard 
state to the so-called high/soft-state (see e.g. McClintock 
& Remillard 2004). Then the X-ray luminosity is dom- 
inated by a strong thermal component originating from 
a geometrically thin, optically thick disc. Current obser- 
vations show that the power-law component is steeper 
(r ~ 2.5 — 3) and much less luminous than in the hard 
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state, suggesting the disappearance of the Comptonising 
plasma in this state. 

Recent multi-wavelength observations of accreting 
black holes in the hard state have shown the pres- 
ence of an ubiquitous flat-spectrum radio emission (see 
e.g Fender 2004), that may extend up to infrared and 
optical wavelengths. The properties the radio emis- 
sion indicate it is likely produced by synchrotron emis- 
sion from relativistic elect rons in compact, self-absorbed 
jets iBlandford fc Konigll 1979; iHiellming fc Johnstori 
1988). This idea was confirmed by the discovery of a 
continuous a nd steady milli arcsecond compact jet around 
Cygnus X-l dStirling et al] 2001V Moreover, in hard state 
sources a tight correlation has been found between the 
hard X-ray and radio luminosities, holding over more 
than three decades in luminosity (Corbel et al. 2003; 
Gallo, Fender & Pooley 2003). In contrast, during high- 
soft state episodes the sources appear to b e radio weak 
(Tananbaum et al. 1972; Fender et al. 1999: ICorbel et al] 
2000), suggesting that the Comptonising medium of the 
low/hard state is closely linked to the continuous ejection 
of matter in the form of a small scale jet. 

Merloni & Fabian (2001) have pointed out that the 
energy content of the electrons in the Comptonising 
medium is too low to account for the observed X-ray lu- 
minosities. The Comptonising electrons have to be tightly 
connected to an energy reservoir where large amount 
of accretion power is stored before being transfered to 
them, and ultimately radiated. Independently, the rapid 
X-ray variability in Cyg X-l also suggests the presence 
of of energy reservoirs (Negoro et al. 1995; Maccarone & 
Coppi 2002). Because angular momentum transport in 
accretion flow is most likely due to magneto-rotational 
instability (MRI, see Balbus & Hawley, 1998), a natu- 
ral candidate for the energy repository is the magnetic 
field amplified in the disc by the MRI-turbulent flow. 
However, if the dissipation of the tangled field (via mag- 
netic reconnection) and non-adiabatic turbulent heating 
preferentially energize the protons, rather than the elec- 
trons, the hot Comptonising plasma could become two- 
temperature, an d the protons themselves act as th e main 
energy reservoir dPi Matteo. Blackman fc Fabia n. 1997). 

Models and simulations of jet production (Blandford 
& Znajek 1977; Blandford & Payne 1982; Meier 2001) in- 
dicate that jets are driven by the poloidal component of 
the magnetic field. Therefore storage of energy into mag- 
netic structures driving the jet and powering the Comp- 
tonising electrons is the most straightforward explanation 
for the observed corona-jet association 1 . In this context, 
the corona would constitute the location where the jet is 
launched. These idea were developed by several authors 
in the context of geometrically thin and/or thick accre- 
tion flows (Meier 2001; Livio, Pringle & King 2003) as 
well as accretion disc coronae (Merloni & Fabian 2002). 

Besides the radio/X-ray correlation observed on long 
(> 1 day) timescales, there are indications that the 
corona-jet coupling operates on timescales as short as 



1 See however MarkofT, Falcke & Fender (2001) and 
Georganopoulos, Aharonian & Kirk (2002) for alternative in- 
terpretations involving dominant X-ray emission from the jet. 



a few seconds or less. The best example is provided by 
the X-ray nova XTE J1118+480 (Remillard et al. 2000; 
McClintock et al. 2001a). During its outburst in 2000, 
this black hole showed all the X-ray properties of hard 
state sources. Fast optical and UV photometry has shown 
rapid optical/UV flickering presenting complex correla- 
tions with the X-ray variability (Kanbach et al. 2001; 
Hynes et al. 2003, hereafter K01 and H03 respectively). 
This correlated variability cannot be caused by reprocess- 
ing of the X-rays in the external parts of the disc. Indeed, 
the optical flickering occurs on average on shorter time- 
scales than the X-ray one (K01), and reprocessing mod- 
els fail to fit the complicated shape of the X-ray/optical 
cross correlation function (H03). Spectrally, the jet emis- 
sion seems to extend at least up to the optical band 
(McClintock et al. 2001b; Chaty et al. 2003, hereafter 
C03), although the external parts of the disc may pro- 
vide an important contribution to the observed flux at 
such wavelengths. The jet activity is thus the most likely 
explanation for the rapid observed optical flickering. For 
this reason, the properties of the optical/X-ray correla- 
tion in XTE J1118+480 might be of primary importance 
for the understanding of the jet-corona coupling and the 
ejection process. 

The simultaneous optical/X-ray observations are de- 
scribed at length in a number of papers (K01; Spruit & 
Kanbach 2001; H03; Malzac et al. 2003, hereafter M03). 
As discussed in these works, the observations are very 
challenging for any accretion model. The most puzzling 
pieces of evidence are the following: (a) The optical/X- 
ray Cross-Correlation Function (CCF) shows the optical 
band lagging the X-ray by 0.5 s, but with a dip 2-5 sec- 
onds in advance of the X-rays (K01); (b) The correlation 
between X-ray and optical light curves appears to have 
timescale-invariant properties: the X-ray/optical CCF 
maintains a similar, but rescaled, shape on timescales 
ranging at least from 0.1 s to few 10 s (M03); (c) The 
correlation does not appear to be triggered by a single 
type of event (dip or flare) in the light curves; instead, 
as was shown by M03, optical and X-ray fluctuations of 
very different shapes, amplitudes and timescales are cor- 
related in a similar way, such that the optical light curve 
is related to the time derivative of the X-ray one. In- 
deed, in the range of timescales where the coherence is 
maximum, the optical/X-ray phase lag are close to 7r/2, 
indicating that the two lightcurves are related trough a 
differential relation. Namely, if the optical variability is 
representative of fluctuations in the jet power output Pj , 
the data suggest that the jet power scales roughly like 
Pj oc — where P x is the X-ray power. 

Here we will show that, if indeed there is a common 
energy reservoir feeding both the jet and the corona, this 
differential relation is naturally satisfied, provided that 
the jet power dominates over the X-ray luminosity. 

We will first present the energy reservoir model and 
suggest a simple physical scenario for the energy reservoir 
and jet disc/coupling (Section |5J. We will then present 
a time dependent model that captures the main fea- 
tures of the multi-wavelength variability observed in XTE 
Jl 118+480, and discuss our main results obtained by 
comparing the properties of the simulated lightcurves 
with the observations (Section [3J • Section 0] will be de- 
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voted to a discussion of the constraints on the nature 
of the accretion flow in XTE J1118+480 derived from 
both spectral and temporal analysis, while in Section 
we make an attempt to generalize these results to other 
accreting black hole sources. Finally, we summarize our 
conclusions in section |S| 



2 THE ENERGY RESERVOIR MODEL 
2.1 A simple analogue 

The time-dependent model that we have developed is 
complicated in operation and behaviour. In order to gain 
some insight into its operation consider a simple model 
consisting of a tall water tank with an input pipe and 
two output pipes, one of which is much smaller than the 
other. The larger output one has a tap on it. The flow in 
the input pipe represents the power injected in the reser- 
voir Pi, that in the small output pipe the X-ray power P x 
and in the large output pipe the jet power Pj . 

If the system is left alone the water level rises until 
the pressure causes P = Pj + P x . Now consider what 
happens when the tap is opened more, causing Pj to rise. 
The water level and pressure (proportional to E) drop 
causing P x to reduce. If the tap is then partly closed, 
the water level rises, Pj decreases and P x increases. The 
rate P x depends upon the past history, or integral of Pj. 
Identifying the optical flux as a marker of Pj and the X- 
ray flux as a marker of P x we obtain the basic behaviour 
seen in XTE J1118+480. 

In the real situation we envisage that the variations 
in the tap are stochastically controlled by a shot noise 
process. There are also stochastically-controlled taps on 
the input and other output pipes as well. The overall be- 
haviour is therefore complex. The model shows however 
that the observed complex behaviour of XTE J1118+480 
can be explained by a relatively simple basic model in- 
volving several energy flows and an energy reservoir. 



2.2 A magnetic energy reservoir ? 

This simple model is largely independent of the physical 
nature of the energy reservoir. In a real accretion flow, the 
reservoir could take the form of either electromagnetic 
energy stored in the X-ray emitting region, or thermal 
(hot protons) or turbulent motions. The material in the 
disc could also constitute a reservoir of gravitational or 
rotational energy behaving as described above. 

In order to be more specific, let us outline the pos- 
sibility of jet disc coupling through a magnetic energy 
reservoir. For the sake of simplicity, we assume that the 
main energy reservoir for the radiating electrons is indeed 
the magnetic field. This would correspond to a system in 
which the heating of the protons by the MHD turbulence 
is negligible. Quataert (1998) and Quataert and Gruzinov 
(1999), have shown that this is the case if magnetic pres- 
sure is not much smaller than gas pressure (low plasma 
f3 parameter), and we will therefore assume this is in- 
deed the case here. Then, the power channelled into the 
particle heating (through magnetic reconnection) and es- 
caping the corona as X-ray radiation can be written as 



P x = (v d /IQE, 



(1) 



where Vd is the field dissipation speed, which depends on 
the details of the dissipation process, P x is the typical 
size of an X-ray emitting region and E = (B 2 /8ty)V is 
the total magnetic energy contained in the hot phase (V 
is its volume). 

The jet power is instead related to the poloidal 
field component. For the purpose of simple esti- 
mates, the MHD jet power can be written as 
llLivio. Qgilvie. fc Pringld 0*99)): 



Pj = (Pp/87r)APA 



(2) 



where A is the area of the disc threaded by the poloidal 
field and R C Q is the typical rotational velocity of the 
field lines. If the field in the disc is amplified by MRI 
turbulence and dissipates mainly in the hot coronal phase 
(i.e. it does not possess a large-scale external component) 
, then its poloidal component can be expressed as B p ~ 
hB, where h = H/R c is the scaleheight of the hot phase 
(Livio et al. 1999; Meier 2001; Merloni & Fabian 2001). 
The jet power can then be rewritten: 



A o 

—h 2 R c nE. 



(3) 



For the sake of simplicity, we assume that the jet 
power is taken mainly from the fangled magnetic field. 
This implies that the outflow is accelerated by magnetic 
dissipation processes and/or that the magnetic energy of 
the corona is carried away by an essentially electromag- 
netic outflow. In fact, it is possible that the jet is powered 
directly by the disc rotational energy without significant 
field dissipation. This would be the case whenever the 
poloidal field is coherent on large enough scales to ex- 
ert a significant torque on the disc (see e.g. the model of 
King et al. 2003). However, such a mechanism for angular 
momentum transport in the disc competes with MRI, so 
that when large scale fields extract energy and angular 
momentum from the disc, the MRI dynamo switches off, 
also draining the magnetic energy reservoir. 

Summarizing, we will make the assumption that 
both the jet and X-ray power are tapped from the same 
magnetic energy reservoir. Note that from equations 
and both jet and X-ray powers scale linearly with the 
total reservoir energy. The total power extracted from the 
magnetic field in the hot phase is then: 



Ah 2 R c Q 



(4) 



If we define the relative fractional power of the jet 
Tj = P\j (P x + Pj), we then find that (see also 
iMerloni fc Fabianl 2002) 

V 



y 7 = 1 



Vd -h- 2 



AR X R C Q 



(•>) 



The discussion so far makes no distinction between 
a corona made of a collection of active phases atop a cold 
disc and a continuous hot inner flow. Indeed the differ- 
ences show up in the geometrical factors in the above 
equations. In the case of a structured corona, we con- 
sider N cylindrical active region (magnetic tubes) located 
atop the disc at distances ranging from 3Ps to P c . The 
typical radius R a and height H a of the active regions 
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scale linearly with their the distance R from the black 
hole, we note a = R a (R)/R and h = H a (R)/R. We 
assume that their radial distribution scales like 1/R so 
that the covering factor of the corona is independent of 
distance. Then the average radius of an active region is 
P x = aP c /ln(P c /3P s ), A = (7v/2)Na 2 R 2 /ln{R c /3R s ), 
and V = (Tv/S)Na 2 hR 3 c /\n{R c /3R s ). 

On the other hand, in the thick disc case, we have 
i? x = R c , A = -kRI and V = (2/3)irhR%, where here R c 
is taken to coincide with the outer radius of the thick disc 
and a constant H/R (wedge- like geometry) is assumed. 



2.3 Time dependent model 

In a stationary flow, the extracted power Pj + P x would 
be perfectly balanced by the power injected into the 
magnetic field, which is, in the most general case, given 
by the difference between the accretion power and the 
power advected into the hole and/or stored in convec- 



tive motions 



Mc z 



Padv ,< 



However, obser- 



vations of strong variability on short time scale clearly 
indicate that the heating and cooling of the X-ray (and 
optical) emitting plasma are highly transient phenomena, 
and the corona is unlikely to be in complete energy bal- 
ance on short timescales. We therefore introduce a time- 
dependent equation governing the evolution of its total 
energy E: 



E = Pi-R- P x , 



(6) 



and we assume that all the three terms on the right hand 
side are time dependent. 

Furthermore, we assume that the optical light comes 
mainly from synchrotron emission in the inner part of the 
jet. A contribution from X-ray light reprocessed in the 
external part of the disc could be present in the optical 
emission. However, this component is likely to be weak 
(cf. KOI; H03). 

The (time averaged) total observed optical flux is: 



O p t oc Pj + / r P x 



(7) 



where the normalization factor / r depends on the geome- 
try of the system and on the optical wavelength, with / r 
increasing towards shorter wavelengths 3 . 



2 In fact, the secular evolution of Convection Dominated Ac- 
cretion Flow leads to accumulation of mass at the outer disc 
boundary, and thus to a non stationary flow. However, if wc 
limit ourselves to observations which are shorter than a typical 
outburst duration, we can regard convection as an alternative 
escape route for the graviational energy dissipated by the ac- 
creting matter out of the region of interest (where X-ray and 
optical emission are produced). 

3 In the fully time dependent model we will consider, addi- 
tional effects on the properties of the optical emission should 
be considered: First of all there is a time delay A between 
matter ejection and dissipation in the form of optical photons. 
As the jet is relativistic and the optical light is produced at 
short distance (< 1000GM/c 2 ) from the hole, the delay should 
be A < 0.1 s. Similarly, reprocessing will also introduce a time 
delay of order of 0.1 s and any variability of P x on time scales 
shorter than a few seconds will not be apparent in the repro- 
cessed light 



We introduce the instantaneous dissipation rates Kj 
and K x : 



P- 3 (t) = Kj(t)E(t) 
P*(t) = K K {t)E(t), 



(8) 
(9) 



and we consistently denote (Kj), (K x ), (Pi), and {E) as 
the time averaged values. We define the average dissipa- 
tion time: 



T dis = ((K j ) + (K x ))- 1 , 



(10) 



such that if the energy reservoir is not fed (i.e. Pi = 0) 
its level decays with an e-folding time Tdi s . Also, it is 
convenient to define 



V = (Kx)T dia 



(11) 



which is the average fraction of the total power that goes 
into the X-ray emission. In the framework of the magnetic 
reservoir of section |2~!H the average values of the dissipa- 
tion rates are related by eqs. Q and to the (time time 
averaged) values of the physical parameter of the system. 
Namely: (Kj) = (Ah 2 /V)R c n and (if x ) = v d /R*. 

Obviously, the detailed physical modeling of the time 
evolution of such a jet corona system is an extremely com- 
plex problem that is beyond the present computer capa- 
bilities. Instead we will adopt a more phenomenological 
approach. We will model the variability of the source by 
assuming random fluctuations of Kj, K x and P, then 
compare the results with the observations and constrain 
the properties of these fluctuations. 

In general, both the instantaneous injected power P 
and the dissipation rates Kj and K^ may depend on the 
amount of energy E stored in the reservoir. However, the 
fact that P and K x are required to vary on time scales 
different from those of P x (see below), suggest that they 
are rather independent of E. Moreover, observations show 
that in black hole binaries and Seyfert galaxies the X-ray 
amplitude of variability is linearly related to the X-ray 
flux level (Uttley & Mc Hardy 2001). In general, if P, 
K x and Kj are either strongly dependent on E, or corre- 
lated with each other, this would introduce a non-linear 
relation between RMS amplitude and observed flux. The 
observations thus suggest that the values of P and the 
dissipation rates Kj and K* at a given time are nearly 
independent of the level of the reservoir energy (although 
E does depend on the history of Kj and K x ). For these 
reasons, as well as for the sake of simplicity, we will then 
assume that Kj, K x and P fluctuate randomly and in- 
dependently, driving the fluctuations of the energy reser- 
voir. 

For the specific form of the fluctuation, we use expo- 
nentially rising shot profiles: 

This profile was chosen for its simplicity. Our results re- 
garding the optical/X-ray correlation are not sensitive 
to the shape of the individual flares. On the other hand 
they will strongly depend on the amplitude, time-scales 
and occurence rate of the shots. The amplitude A and 
the occurrence rate A of the random shots are taken con- 
stant. These quantities are related to the average dissi- 
pation rate. The shot duration r is distributed within 
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Train and T ma x with a power law distribution p(r) oc t~ p . 
These parameters constrain the fractional amplitude of 
variability and the shape of the power spectrum of the 
fluctuations that are imposed on the system. The frac- 
tional RMS scales like 1/VA, and the power spectrum is 
a power law ranging from frequencies 1/T max to l/r min 
with a slope a — 3 — p (see e.g. Poutanen and Fabian 
1999). For K x , K, and Pi, the parameters A, p, T m i n , 
Tmax are, in general, not identical. In presenting our re- 
sults, the subscripts x, j and i will refer to A" x , Kj and 
Pi respectively. 

For a specific set of parameters we first generate time 
series for K x , Kj and P, solve the time evolution of the 
energy reservoir E and then use it to derive the the result- 
ing optical and X-ray light curves. The ejected material 
travels from the corona to the shock region or photo- 
sphere where the optical photons are produced. As dis- 
cussed above, this introduces a time delay of A < 0.1 s in 
the optical emission. To model this delay, we simply shift 
the optical light curve by A. We then compute the X-ray 
and optical Power Density Spectra (PDS), Auto Correla- 
tion Functions (ACF), their Cross-Correlation Function 
(CCF), coherence and phase lag spectrum for compari- 
son with the observed ones, in order to obtain the largest 
possible number of different observational tests for our 
variability model. 

2.4 Main observational constraints 

Before we proceed and illustrate the results of our simu- 
lations, a remark is in place with respect to the require- 
ment any specific realization of our model will have to 
fulfill in order to reproduce the main observational char- 
acteristics. This will guide us in the exploration of the 
vast parameter space of the model, as well as provide us 
with useful insight on the physical interpretation of the 
results. 

We first note that combining equations © and © 
we obtain for the total instantaneous jet power the fol- 
lowing relation: 

P:=Pi-(l + J)Px-Px/K x . (13) 

We can see from this equation that the differential 
scaling Pj oc -P x , observed in XTEJ1118+480, will be 
rigorously reproduced provided that: 

• i^x is a constant; 

• Pi — P x is a constant. 4 

It is physically unlikely that those conditions will 
be exactly verified. In particular, P x is observed to have 
a large RMS amplitude of variability of about 30 per- 
cent. However, the observed differential relation holds 
only roughly and only for fluctuations within a relatively 
narrow range of time-scales 1 — 10s. Therefore, the above 
conditions need only to be fulfilled approximatively and 
for low frequency fluctuations (> Is). In practice, the 

4 This constant can differ from because the differential scal- 
ing is observationally demonstrated only for the varying frac- 
tion of the optical and X-ray fluxes. 



following requirements will be enough to make sure that 
the low frequency fluctuations of the right hand side of 
equation 1131 are dominated by P x : 

• P x <C Pi, implying that the jet power, on average, 
dominates over the X-ray luminosity; 

• the amplitude of variability of K x and P in the 1-10 
s range is low compared to that of Pj . In other words the 
1-10 s fluctuations of the system are mainly driven by the 
jet activity. 

The first condition is crucial, and will be discussed 
in more detail in Section^] The second condition requires 
that the mechanisms for energy reservoir filling and dissi- 
pation in the corona and in the jet are occurring on quite 
different timescales, and this will also provide us with ad- 
ditional constraints on the global dynamical properties of 
the system. 



3 RESULTS 

In figure we show the results of a simulation with 
Pdis = 0.5 s and / x = 0.1 (see TableQfor the value of the 
other parameters). The model produces an X-ray power 
spectrum with a plateau up to ~ 0.1 Hz and a power- 
law component with slope ~ 1.4 above that frequency, 
with most of the X-ray variability occurring around 0.1 
Hz. The optical PDS power-law has a flatter slope (~ 1) 
up to 1 Hz and then softens to a slope similar to that 
of the X-ray PDS. The resulting optical ACF is signif- 
icantly narrower than the X-ray one. The full-width-at- 
half-maximum (FWHM) of the two ACFs differs by a 
factor > 2. The overall coherence is low (< 0.4): reaching 
a maximum in the 0.1-1 Hz range and decreasing rapidly 
both at lower and higher frequency. The phase-lags are 
close to 7r/2 in the 0.1-1 Hz range and increase from 
at low frequencies up to ir at around 6 Hz. At higher fre- 
quencies the phase lags spectrum is characterized by large 
oscillations. Finally the resulting CCF rises very quickly 
at positive optical lags, peaks around 0.5 s (this is the 
post-peak) and then declines slowly at larger lags. The 
two bands appear to be anti-correlated at negative opti- 
cal lags indicating a systematic optical dip 1-2 s before 
the X-rays reach their maximum (pre-dip). 

All these characteristics are observed in XTE 
J1118+480. Moreover, the integrated X-ray RMS is 25 
percent as observed. The optical variability is 20 percent 
slightly larger than the 16 percent reported by M03 but 
this could be reduced if the expected constant disc com- 
ponent was added. 

Obviously the model parameter values were carefully 
selected in order to reproduce these characteristics. Given 
the complexity of the model, a full investigation of the 
parameter space is premature. However, it may be useful 
to illustrate the main effects of the different parameters 
by considering several simple situations. 
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3.1 Variability dominated by the jet dissipation 
rate Kj 

3.1.1 Jet dominated models (fx«l) 

We first consider the case where the X-ray emission is 
energetically negligible. The jet fully drives the variabil- 
ity of the system in the limit of large Ai and A x , and 
low / x . In this limit, both conditions for the applicability 
of the model discussed in the previous section are ex- 
actly verified. We thus expect the scaling Pj oc — P x to be 
realised. Indeed, we obtain phase lags that are close to 
7r/2 independent of frequency. In this limit, the resulting 
CCF is purely anti-symetric with pre-dip and post-peak 
of identical amplitude and lag (contrary to what is ob- 
served). The optical lag of the post-peak is controled by 
Tdi s , which sets the time scale on which the energy reser- 
voir responds to the imposed fluctuations. At frequencies 
larger than 1/Tdi s the X-ray variability decreases because 
there cannot be any variability of the reservoir energy 
E on time scales larger than Tdi s . The coherence func- 
tion also decreases at frequencies above 1/Tdi s . On the 
other hand, at frequencies lower than 1/Tdi s the optical 
variability decreases because at such low frequencies the 
system evolves in a quasi-static way. Since Pi is constant 
and P x is energetically negligible, the variability of Kj is 
compensated by the variability of E so that Pj remains 
almost constant, too. 



3.1.2 X-ray dominated models (fx~l) 

If / x is large but the amplitude of variability of K x and 
Pi are kept negligible by increasing A x and Ai, then the 
power output can be dominated by the X-rays but the 
variability is still driven by the jet. 

In this case, we still have n/2 phase-lags but only 
at frequencies > 1/Tdi s . The effects of a large X-ray dis- 
sipation is the appearance of phase-lags -k at frequen- 
cies < 1/Tdis indicating an anti-correlation. This anti- 
correlation is also apparent in the CCF, which tends to 
be negative and symmetric around zero lag. The anti- 
correlation at low frequencies is due to the fact that on 
long time-scales the system is always in equilibrium, but 
contrary to the jet dominated case, X-ray losses are not 
negligible in this case. The reservoir energy E changes 
in order to keep the total output constant (and equal 
to P), but Pj does change. For this reason the fluctua- 
tions of the jet output are strictly anti-correlated to that 
of the X-ray power. Such an anti-correlation is in clear 
conflict with the data, the observed phase-lags indicating 
rather a correlation at low frequencies. Thus, even if the 
jet drives the variability as in the previous example, we 
definitely need / x <C 1 in order to reproduce the data. 
In practice, exploring a large volume of the parameter 
space we could not find a reasonable agreement with the 
data for / x larger than ~ 0.2. This is a very robust result 
of our time-dependent modeling, and we will examine its 
consequences in Section 2] 



3.2 Variability dominated by the coronal 
dissipation rate K^ and/or by the power 
input Pi 

The situation where the variability is dominated by the 
X-ray dissipation rate in the corona is perfectly symmet- 
ric to the case where the variability is dominated by the 
jet, studied above, with however an inversion of the sign 
of phase lags and time-axis in the CCF plot. 

On the other hand, fluctuations of P introduce per- 
fectly positively correlated fluctuations of Pj and P x . If 
the fluctuations of P dominate we obtain X-ray and op- 
tical light curves that are perfectly correlated on all time 
scales with identical power spectrum and unity coher- 
ence, in contrast with the observations. 

In general, when in a given range of timescales more 
than one component (i.e. P, Kj or K*) drives the vari- 
ability, the main effect is a strong reduction of the coher- 
ence function in the corresponding range of frequencies. 

3.3 A realistic case 

Now that we understand the basic effects of the different 
parameters, we get back to our 'realistic' model and ex- 
plain how we were guided toward this solution and how it 
reproduces qualitatively all the timing features observed 
in XTE Jl 118+480. 

As discussed above, reproducing the data requires 
the jet to dominate the energy budget and be the main 
driver of the variability (at least for the 1-10 s fluctua- 
tions). Therefore we set / x =0.1 and the parameters Xj, 
Pj, T m in j and T ma x j were adjusted in order to have large 
enough 1-10 s fluctuations of the jet power. These pa- 
rameters also control the shape of the optical power spec- 
trum. In particular, the slope of the shot distribution pj 
is chosen so that it leads to the observed slope of the op- 
tical power spectrum above 1 Hz. The dissipation time is 
set to Tdi s = 0.5 in order to match the observed optical 
lag in the CCF reported by K01. 

The values of Ai, pi, T mini and Tm aX i were fixed so 
that the variability of P introduces some correlated vari- 
ability at low frequencies, enabling us to reproduce the 
asymmetry of the CCF and the small phase-lags observed 
at low frequency. The loss of coherence below 10 Hz is due 
to the simultaneous large variability of P and Kj at low 
frequencies. 

Finally, we also introduce fluctuations of the X-ray 
dissipation rate at high frequencies. This enables us to 
reproduce the loss of coherence observed above one Hz. 
Also X x , p x , T min x and T m ax x were chosen in order to 
reproduce the X-ray power spectrum at high frequencies 
which is otherwise too steep. 

A propagation lag of A =0.05 s then provides the 
observed high frequency behaviour of the phase lags. If 
the propagation lag is neglected, the phase lag spectrum 
is flat at high frequencies with phase lag at a constant 
value ~ 7r/2. The constant time delay introduces an ad- 
ditional phase lag ip = 2nfA. As the time-scale of the 
fluctuations approaches that of the time delay this ad- 
ditional phase lag becomes dominant. The overall phase 
lag thus increases and reaches n at ~ 6 Hz. The phase 
lags are defined only in the range [— tv,tv] and can be 
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only measured modulo 2n. Thus, at higher frequencies 
the phase lag shifts to -tt then increases linearly to n 
and so on, producing large 'oscillations' in the pahse-lag 
spectrum. A similar behaviour of the phase lags above 
1 Hz is observed in XTE J1118+480 (see M03, Fig. 5), 
indeed suggesting a propagation lag of A ~ 0.05 s. On 
the other hand, the overall effects of propagation delays 
on the CCF and coherence function are negligible. 



3.4 The effects of reprocessing 

We now consider the effects of a possible reprocessing 
component from the disc. We model this component by 
convolving the X-ray lightcurve with a transfer function 
T r describing the time delay and variability smearing due 
to reprocessing. We use the following transfer function : 



T r (t) 



oc - exp(— t) exp 
o 



+ 0.28 



exp 



t 

~ 5 

ex P r-(V5) 20 l (14) 



where t is measured in seconds and the reprocessing time 
delay 5 = 0.1 s. 

This function is a rough analytic approximation to 
the theoretical transfer functions shown in H03. 

The reprocessed component is added to the optical 
light curve. Its normalization is parametrized by the ratio 
of the average flux of the reprocessed component to that 
of the jet component. In the simulation shown in Figure|5| 
we take this ratio to be 0.5, keeping all other parameters' 
values identical to those of the previous simulation. The 
main effect is to add a low frequency component to the 
optical lightcurves highly correlated to the X-ray one. By 
comparison with Figure the optical ACF is broader, 
the coherence is higher at low frequency and the optical 
pre-dip tends to disappear, the optical lags are shorter in 
particular at low frequency. 

The increase of the reprocessed component is asso- 
ciated with an increase in the coherence function, and 
the X-ray and optical ACF and PDS becoming similar 
and shorter optical lags. We expect to observe this evolu- 
tion if we look at shorter and shorter wavelengths where 
the disc emission increasingly dominates the total op- 
tical/UV flux. Indeed H03 report similar trends in the 
dependence of the correlation upon wavelength. We also 
note that at low frequencies (< 0.1 Hz), H03 observe a 
large coherence in the far UV domain, while M03 report 
almost zero coherence. This could be interpreted as the 
effects of reprocessing in the far UV that would become 
negligible at optical wavelengths. 



4 THE NATURE OF THE 

LOW-LUMINOSITY ACCRETION FLOW 
IN XTE J1118+480 

4.1 Constraints from the time variability study 

The detailed analysis of the complicated structure of the 
multi-wavelength variability of the source have allowed 
us to select a region of parameter space for which our 
model successfully reproduces all the main observational 



-Mis — u.o s 


f — n i 
jx — U.J 




Aj = 100 s^ 1 


Aj = 50 s" 1 


A x = 1000 s _1 


Pi = 2.1 


Pj = 1-4 


Px = 1.1 


^mm i — 0.1 S 


T~min j — 0.01 S 


^mm x — 0-01 S 


Tmax i — 7 S 


7~max j — 10 S 


1~max x — 0.5 S 


RMS P ; = 0.20 


Kj = 0.27 


Kx = 0.11 


RMS E = 0.23 


Opt = 0.21 


X = 0.26 



Table 1. The model parameters used in the simulations of 
Fig. The two last rows show the resulting fractional am- 
plitude of the dissipation rates, energy reservoir, optical and 
X-ray fluxes. In model|2]some reprocessing was added and the 
optical fractional RMS is then reduced to 0.16. 



pieces of evidence. The parameters are summarized in 
Table Q In principle these parameters can contrain the 
dynamic and geometry of the accretion flow. However this 
would require a detailed physical model for the jet/disc 
coupling. Here we will assume that, as suggested in sec- 
tion the magnetic field constitutes the main energy 
reservoir and investigate the consequences for the accre- 
tion flow. 

Let us start by examining the two parameters that 
reflect the time-averaged global energetics of the system: 
Tdis and / x . From Eqs. Q, @, l(Tfljl and ifTT) . we see that 
they are related to the physical parameters of the system: 

Vd _ fx 

Rx Tdis 



(15) 



— ft vk{Rc 



fx 



Tdis 



(16) 



where we have assumed that the disc rotation law is Ke- 
plerian and we denote Vk(R c ) the Keplerian speed at the 
radius R c . 

In the estimates below, we will assume, as we do 
throughout the paper, a black hole mass of WMq for 
XTE J1118+480. Equation QSJ is a constraint on the 
magnetic field dissipation speed in the hot phase: for f K = 
0.1 and Tdis = 0.5 it can be rewritten as: 

-^2x 10- 5 r x , (17) 
c 

where r x = R x /Rs- 

On the other hand, Equation 1161 constrains the 
corona/hot flow scaleheight, which should be defined by 
selecting one of the two possibility for the hot phase ge- 
ometry. It turns out that for both geometries we obtain 
the same relation: 



ft ~ 1.7 x 10" 



0-5 y f r c \ 



3/2 



(18) 



T dis 0.9 U00/ 

where r c = R c /Rs- 

In the case of the thick inner flow, self-consistency 
would require a substantial scaleheight of the hot flow. 
If, for example, we require ft > 0.1, we obtain a con- 
straint on the radial extent of the inner hot flow r c > 70, 
which, incidentally, appear to be consistent with the 
spectroscopically-inferred inner radius of the truncated 
cold disc (Esin et al. 2001; Chaty et al. 2003). 

The same spectroscopic constraints also apply to the 
patchy corona case (see next section) and the radial ex- 
tension of the corona above the disc should also be at 



© 2002 RAS, MNRAS 000, [I]-?? 



8 Malzac, Merloni & Fabian 



1.5 
0.5 

1.5 

K? 1 
0.5 

1.5 

k H 1 
0.5 



1.5 
H 1 
0.5 

1.5 

5- 1 

0.5 





(a) 





(c) 



100 



102 



104 106 
time (s) 



108 



110 



1.0 

C 0.8 

o 
+j 

S3 

t 0.6 
u 
u 
o 

g 0.4 

3 

" 0.2 
0.0 
0.3 

J 

I 0.2 

J 

, o.i 



0.0 



a. 
o 



-0.1 



_ — — - optical 




! X-rays 

'- / 




: J i 
- / ' 


i\ : 
i \ j 
\ \ 

\ \ : 
\ \ 













s 

CVJ 

c 



O 

a, 





Pi 

Kj 

K x . 


■ i' / / "' 




| (b) :'V 


\ s 




-10 



10 



Opt. lag, (s) 



-1 1 

log Frequency, Hz 



Figure 1. Results of simulations for parameters given in tablc^ Sample time series (panel a) and power spectra (panel b) of 
Pi, Kj, K x , resulting X-ray and optical fluxes light curves (panel c), X-ray /optical autocorrelation and cross-correlation functions 
(panel d), power spectra, coherence and phase-lags (panel e). 
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Figure 2. Same simulation as in Fig. Q but adding a disc reprocessing component accounting for 30 percent of the optical flux 
on average (see section [3.41 . Panel a: sample light curves of the reprocessed component and total optical Panel b: X-ray /optical 
autocorrelation and cross-correlation functions. Panel c: power spectra of the reprocessed component, total optical and X-rays 
(top), X-ray/optical coherence and phase-lags (bottom). 



least r c ~ 100. Then for typical inferred values of the 
size of the active regions (r x ~ a few, see e.g Haardt, 
Maraschi & Ghisellini 1994; Di Matteo, Celotti & Fabian 
1999), the aspect ratio of the cylindrical regions writes: 



^1.6 



_05__[L-3- /^\ 5/2 ln(r c /3) 
T dis 0.9 r x U00/ ln(100/3)' 



(19) 



It appears to be slightly larger than unity, as required to 
reproduce the hard X-ray spectrum (see e.g. Malzac, Be- 
loborodov & Poutanen 2001). Note that the aspect ratio 
is very sensitive to the coronal radius. If r c is as large as 
350, as inferred by Chaty et al. (2003) we obtain a quite 
large ratio h/a ~ 27, then internally generated cyclo- 
synchrotron radiation should dominate over reprocessing 
as a source of seed photons for the Comptonisation pro- 
cess. 

Moreover, combining the expression that relates the 
observed X-ray luminosity (0 with the constraints ob- 



tained above on the field dissipation speed, we can esti- 
mate the value of the magnetic field. In the case of the 
thick inner flow we obtain: 



B ~ 3.2 x 10' 



~ 1/2 Tdte f r c \ 
0.5s VIOOJ 



9/4 



(20) 



0.1 0.9. 

From the lower limit on the inner disc radius derived 
above, we conclude that B < 7 x 10 6 G. In the case of 
the patchy corona: 

1 -1/2 

2.5 x 10 7 



B 



/ x n ln(r c /3) N 



0.1 0.91n(100/3) 10 



5/4 



(21) 



T dis 3 ( r c \ 
0.5s r x UOOy 
We note that the RMS variability amplitude should 
scale as iV~ 1//2 , so that the magnetic field intensity is 
directly proportional to the observed variability level. 
Finally, let us discuss the different variability 



© 2002 RAS, MNRAS 000, [I]-?? 



10 Malzac, Merloni & Fabian 



timescales of the dissipation rates, as fixed in our fidu- 
cial model. On the one hand, the range of timescales on 
which the variability of is maximal are between 50 
and 1 Hz, corresponding to the Keplerian frequencies be- 
tween 10 and 100 R s (for a 10 M e black hole). This 
is consistent with the X-rays being produced in the hot 
phase of the accretion flow (either in a geometrically thick 
disc or in the corona). In general the dissipation of the 
tangled magnetic field into particle heating is essentially 
a local process acting on timescales comparable to the 
local dynamical time. Therefore the minimum and maxi- 
mum variability timescales of if x reflect the smallest and 
largest radii where we have substantial energy in the hot 
phase/corona. 

On the other hand, we see from Figure Q that the 
variability of Kj is high on much wider range of timescales 
than K x , from 0.01 to 10 seconds. As discussed in sec- 
tions 1231 an d ED the data require that at low frequencies 
the fluctuations of Kj dominate over the other sources 
of variability. The physical interpretation is that gen- 
crating a powerful outflow requires large-scale coherent 
poloidal field structures threading a significant fraction 
of the inner disc surface and extending also in the ver- 
tical direction. Such structures may take more time to 
build and destroy. So, locally, the timescale of variation 
of Kj can go from the dynamical time itself to larger 
timescales (in the rare cases where the field builds up its 
poloidal component to a larger scale). Based on simple 
arguments, Livio, Pringle & King (2003) estimate that 
the time-scale for establishing a change in the scale of 
the poloidal component of the magnetic field through dy- 
namo processes is larger than the dynamical time scale 
by a factor tj e t/tdyn ~ 2 1 ' fc . Since our model indicates 
tjet/tdyn ~ 10, this would imply h ~ 0.3, which appears 
to be in agreement with the independent estimate pro- 
vided by equation 1181 

We also note that the variability of Kj induces large 
amplitude fluctuations of the energy reservoir that in 
turn lead to an important X-ray variability on time-scales 
much larger than the dynamical time. In fact, most of 
the observed X-ray variability is caused by the jet activ- 
ity, providing an explanation for the long observed time- 
scales. Recently King et al. (2003) have proposed a de- 
tailed physical model for the X-ray variability of accreting 
black holes that considers a similar situation where large- 
scale ejection events modulate the strong X-ray emission 
from the inner disc. 



4.2 Constraints from time averaged 
spectroscopy 

We have shown in section |3] that there is an important 
condition for our a model to work, namely that the jet 
power should dominate over the X-ray luminosity. We 
will now discuss how realistic this assumption is. 

First, let us examine the observed energetic out- 
put of XTE J1118+480. We stress that this is the best 
source available for such kind of study, given its unsur- 
passed spectral coverage. In the following we will express 
all the luminosities in units of the Eddington luminosity 
(L E dd = 1.3 x 1O 38 (M/M ) erg s"\ with M = 10M Q ). 



The numbers given below are taken from the results of the 
multi- wavelength spectral analysis of Chaty et al. (2003) . 
These authors decomposed the overall SED into three 
components: 

• The jet component has a total observed radiative lu- 
minosity of Lj = 3.8 x 10 -5 , assuming isotropy and in- 
tegrating the power-law emission from the radio to the 
optical band. In the optical/UV the emission should start 
to be dominated by the cold accretion disc, but the pos- 
sibility that the jet spectrum extends up into the UV 
domain cannot be ruled out. Indeed, the fact that H03 
observe similar correlated variability in the UV and opti- 
cal lend support to this hypothesis. If this is the case, the 
estimated jet luminosity can increase by up to a factor of 
8. 

• For the standard cold accretion disc we have: Ld = 
1.7 x 10 -3 . This estimate was obtained by Chaty et al. 
(2003) with a fit of the optical to EUV spectrum with 
a multicolour blackbody model. The inner disc temper- 
ature was found to be very low kT ~ 20 eV, suggesting 
that either the inner disc is truncated at ~ 350 -Rs, or 
that the bulk of the accretion power is driven away in a 
coronal outflow and the disc is left very cold. 

• Finally, for the luminosity of the X-ray emitting 
plasma, or corona, we have: L x = 1.2 x 10~ 3 , assuming 
isotropy. 

Summing up all the three above spectral compo- 
nents, we find that the bolometric luminosity of the 
source is: 



Lboi = Lj + L d + L x ~ 3 x 10" 



(22) 



The total mass accretion rate onto the black hole 
can be estimated from the observed luminosity of the 
cold disc component. We define the Eddington-scaled 
accretion rate m = eMc 2 / L^dd where e is the New- 
tonian accretion efficiency e — 1/12. If the cold, geo- 
metrically thin and optically thick disc is indeed trun- 
cated at Ri n = 350 -Rs, as the spectral analysis of Chaty 
et al. (2003) suggests, then standard accretion formulae 
(Frank, King and Raine 2002) give: 

(23) 



0.2. 



m = 4e— —L 

Rs 

On the other hand, if the accretion disc extends down 
to the innermost stable orbit (Ri n = 3Rs), but is sand- 
wiched by a powerful corona, where a fraction /h of the 
power is dissipated, then we can estimate that 



M = i + 



L d (l-V) 



(24) 



The 7] parameter is only poorly constrained, but in any 
case /h should be very large in order to explain the low 
inner disc temperature (Merloni, Di Matteo & Fabian 
2000). Taking rj ~ 0.9 as indicated by the time depen- 
dent model, we obtain /h = 0.88, and, for the total mass 
accretion rate, 

Rin Ld 



m = 4e 



Rs (1 - /h) 



0.015. 



(25) 



In both cases, the estimated accretion rate is much larger 
than the observed bolometric luminosity. Thus, most of 
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the accretion power is not radiated and the source ap- 
pears to be radiatively inefficient. 

The important issue, however, would be to determine 
whether the missing accretion power escapes the system 
in the (low radiative efficiency) jet or in other forms of 
non-radiative losses, such as a slow wind, or large scale 
convective motions, or advection into the black hole. The 
answer to this question resides in the exact determination 
of the jet kinetic power. Unfortunately, there are major 
uncertainties in this determination, mainly because the 
jet radiative efficiency is not known. The jet is expected 
to be a poor radiator because most of the energy is lost 
in adiabatic expansion, thus, although the radiation from 
the jet represents a small fraction of the bolometric lumi- 
nosity the jet could dominate the energetics. Both theory 
and observations indicate that the efficiency can be ex- 
tremely low (as low as e s ~ 10~ 4 see Celotti and Fabian 
1991), but all present estimates are model dependent. In 
general they indicate radiative efficiencies of the order of 
ej ~ 0.01. For the case of XTE J1118+480 this would 
already imply that the total jet power dominates over 
the X-ray luminosity. Obviously, if the jet efficiency is 
lower, or the jet spectrum extends to shorter wavelength, 
or beaming effects are important, the jet dominance can 
easily be much larger. As discussed above, the analysis of 
our time dependent modeling strongly requires f K < 0.1, 
corresponding to a jet efficiency tj < 3 x 10~ 3 . 

In the case of the truncated disc plus inner hot inner 
flow model, the derived m sets an upper limit to the jet 
power, in the case where all the accretion power is lost in 
the jet: P;/L^ ~ m/L x ~ 200. On the other hand, if the 
jet power is comparatively modest e.g. Pj/L x ~ 10 as as- 
sumed in the simulation of Fig.0 this would imply that, 
although the jet dominates over the X-ray emission, it 
does not dominate as a power sink. If the dominant power 
sink is, for example, advection onto the black hole, the 
mass accretion rate we derive appears somewhat large to 
be consistent with ADAF models (Narayan & Yi 1994). 
However, it is important to notice that the observed ac- 
cretion rate is very sensitive to the value of the inner disc 
radius (see eg . 12311 . which is actually poorly constrained. 
Moreover, the critical accretion rate for ADAF models 
above which the ADAF solution breaks down scales as 
a 2 (Rees et al. 1982; Narayan & Yi 1995), and the un- 
known standard a parameter could be larger than usually 
considered leading to the existence of solutions at higher 
rh than previously thought. Finally, we note that the dy- 
namics of the ADAF solutions coupled with jets have 
not been worked out yet. The structure of the accretion 
flow as well as the (in)stability of the solution at high rh 
are likely to be strongly affected by the presence of the 
jet. Although an advection-dominated accretion flow is 
in many aspects consistent with our analysis, it should 
certainly differ considerably from the standard solution. 

In the case of the accretion disc plus patchy corona 
(or coronal outflow) model the rh derived assuming n = 
0.9 is consistent with all of the non-radiated power driv- 
ing the jet. However, when we attempted to fit the 
optical-UV spectrum of XTEJ1118+480 with the disc 
corona solution of Merloni & Fabian (2002), we found 
that it was not possible to produce inner disc tempera- 
tures as low as observed (~ 20 eV). The coronal outflow 



model could be reconciled with the EUV spectral data 
only if an even larger fraction of the power was extracted 
from the disc into the corona/jet, requiring /h to be very 
close to 1. However, in the framework of magnetic coro- 
nae powered by MHD turbulent disc, having / ~ 1 also 
requires a very large (magnetic) viscosity (Merloni 2003), 
and implies an even stronger jet dominance of the total 
energy budget of the source. 

An alternative way to obtain a low inner disc temper- 
ature together with reasonable value of /h would be that 
the cold disc does not extend down to the last stable or- 
bit. Instead, below a few tens of Rs the disc is truncated 
and the flow becomes very inefficient. However, such a 
picture would be very similar to the truncated disc plus 
hot inner flow scenario discussed above, except that the 
jet would be launched from the external accretion disc 
corona. 



5 SOME WIDER IMPLICATIONS OF THE 
MODEL 

5.1 Comparison with other sources 

The behaviour of XTE Jl 118+480 demonstrates that the 
accretion flow onto a black hole can lead to a strong jet 
with little disk emission. A similar situation occurs in 
the massive elliptical galaxy M87 where powerful jets are 
associated with an otherwise weak galactic nucleus. The 
Bondi accretion rate onto the supermassive black hole in 
M87, deduced from the properties of the surrounding gas 
with Chandra observations (Di Matteo et al 2002), can 
give the observed jet power, deduced from the cavities 
made by them in the surrounding intracluster medium, 
provided that accretion efficiency is 0.1 or more. Assum- 
ing that the flow is in a steady state, such a high ef- 
ficiency argues that most of the accreting gas flows to 
within a few gravitational radii of the black hole. Since 
thick hot flows are unstable to mass loss over a wide range 
af radii ( Begelman & Blandford 1999; Quataert & Gruzi- 
nov 2000; Stone, Pringle Begelman 2000), this favours a 
magnetically-dominated cold disk flow in the case of M87. 
By similarity it argues for a magnetically-dominated cold 
disc in XTE J1118+480 as well. The reservoir and pre- 
sumably the base of the jet must also lie within a few 
gravitational radii of the black hole in both cases. 

As mentioned in the Introduction, the low state of 
many Galactic Black Hole systems appears to contain a 
powerful jet. When systems drop below about one per 
cent of the Eddington accretion rate the disc may be 
magnetically dominated (Merloni & Fabian 2002; Livio, 
Pringle & King 2003), with a fast jet taking most of 
the gravitational energy released. This may also apply 
to massive black holes, as illustrated above with M87. 
The issue of how apparently well-fed massive black holes 
in galactic nuclei are feeble sources of electromagnetic 
radiation has been a puzzle for some years (Fabian & 
Canizares 1988; Fabian & Rees 1995; Di Matteo et al 
2001, 2003; Loewenstein et al 2001; Pellegrini et al 2003). 
Most of these objects do however have jetted radio emis- 
sion (Franceschini, Vercellone & Fabian 1998) and a so- 
lution in which the accretion power is principally car- 
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ried away by jets is a strong possibility. Unless the en- 
ergy contained in black hole spin is a key factor, then 
magnetically-dominated discs should be part of that so- 
lution. 

Further time-dependent studies of the multiwave- 
length behaviour in other stellar-mass and massive 
black hole systems may show similar behaviour to that 
of XTE 1118+480. Some examples which already of- 
fer tantalizing behaviour are the Galactic microquasar 
GRS1915+105 (Mirabel et al, 1998; Klein-Wolt et al., 
2002), the Galactic binary GX339-4 which has also dis- 
played rapid optical flickering (Motch et al 1982) and the 
blazar 3C120 which shows X-ray dips before radio jet 
events are seen (Marscher et al 2002). 



5.2 Jet-dominated sources in the low/hard 
state 

Whatever the actual structure of the accretion flow, we 
have shown that during the outburst of XTE J1118+480, 
the total kinetic jet power should dominate over the 
X-ray luminosity, and could possibly be the dominant 
repository of the accretion power. There are additional 
independent arguments in favour of jet dominance in 
low/hard state sources and in XTE J1118+480 in par- 
ticular. Based on the observed radio flux (Lb.) and X- 
ray correlation observed in hard states sources (Falcke 
& Biermann 1996; Gallo, Fender & Pooley, 2003), as 
well as on standard synchrotron formulae (Heinz & Sun- 
yaev, 2003), Fender, Gallo & Jonker (2003) have shown 
that, provided that advection into the black hole horizon 
and/or convective motions do not store a large fraction 
of the accretion power, there should exist a critical ac- 
cretion rate, m cr , below which an accreting black hole is 
jet-dominated. 

The exact value for the critical accretion rate could 
be inferred from the observations, if we knew the total 
jet power at a certain X-ray luminosity, and is given by 
m cr = 2Pj 2 /L x , corresponding to a critical X-ray lumi- 
nosity L x , cr = m cr /2. Fender et al. (2001) derived a lower 
limit for the jet to X-ray power ratio in XTE J1118+480: 
Pj I L x = 0.2, and FGJ03 used this conservative estimates 
to determine the value of the critical rate m cr ~ 7x 10~ 5 . 
However such a low value of the critical luminosity leads 
to several problems. 

First, as shown in FGJ03, during the transition from 
a disc to a jet-dominated state, the dependence of the X- 
ray luminosity on the accretion rate changes from being 
L x oc rh 2 , the right scaling for radiatively inefficient flows, 
to L x oc rh, the scaling for radiatively efficient flows (see 
Fig. 1 of FGJ03). This would imply that with L x ~ 10" 3 , 
XTE J1118+480 should be a radiatively efficient system. 
As discussed above, there is however strong observational 
evidence of the contrary. 

Second, during the transition from jet to disc dom- 
inated state, the jet power changes from Pj oc rh which 
is the natural scaling expected in most jet models, to the 
seemingly unphysical Pj oc rh 1 / 2 . 

Furthermore, black holes in the hard state should 
show some kind of spectral transition in the X-ray band 
at the critical luminosity Z/ XjCr ~ 3 x 10 -5 , due to the 



drastic changes in emission mechanisms that are needed 
to account for the different scalings of L x with the accre- 
tion rate. The observations of low/hard state sources at 
such low luminosities are few and hard to perform, how- 
ever no indication of any dramatic spectral change in any 
hard state source down to quiescent level has ever been 
reported (Kong et al., 2002; Hameury et al., 2003) 

In fact, the only physical transition that we do ac- 
tually observe is the transition between the hard and the 
soft state that occurs at luminosities of at least a few 
percent of Eddington luminosity (Maccarone 2003). We 
believe that, if the above mentioned difficulties are to be 
solved, then m cr has to correspond to luminosities that 
are comparable to, or larger than, hard-to-soft state tran- 
sition luminosities. 

For the case of XTE J1118+480, instead of using the 
lower limit for the jet to X-ray power (Pj/L x = 0.2), we 
can adopt the much larger value Pj/L x ~ 10 required by 
our variability model. Then we find m cr ~ 0.2, involving 
a transition at L x , cl - = m cr /2 ~ 0.1. This is in agreement 
with the idea that the transition from jet dominated to 
X-ray dominated states occurs at luminosities similar or 
slightly higher than the hard to soft state transition. We 
stress that the arguments of FGJ03 holds only as long as 
the Lr — Z/ x correlation is valid i.e. not in the soft-state. 
There is thus no reason to associate the disc dominated 
sate of FGJ03 to the soft state. 

If, as we suggest, the critical luminosity L XjCr is at or 
above the hard to soft state transition luminosity, then 
the disc dominated state of FGJ03 does not exist. More- 
over, although one can speculate that the transition to 
soft state may be be triggered by the loss of the jet dom- 
inance, there is a priori no reason for the critical lumi- 
nosity L XjC r to be coincident with the hard to soft state 
transition luminosity. In other words, the observed pres- 
ence of a hard-to-soft state transition does not provide 
an upper limit to the jet power, but only a lower limit. 

We also note that the above discussion was made 
neglecting any contribution from advection and other 
possible non-radiative losses. However, in jet-dominated 
sources the jet power scales proportionally to the accre- 
tion rate, as is the case for any adiabatic loss. Therefore 
the effects of advection (and/or convection) would not 
change the scalings of jet and X-ray powers provided that 
we replace the accretion rate rh with rh e s = rh(l — fad), 
where / a( j is the fraction of accretion power that is ad- 
vected into the black hole. In particular, the critical X-ray 
luminosity L x , C r, although corresponding to a larger 'real' 
rh, would remain unchanged. 

To conclude, the central requirement of our model 
i.e. that in XTE J1118+480 the jet power dominates over 
the X-ray luminosity, appears to be supported by the ob- 
servations of hard state black holes. Then, if the argu- 
ments of FGJ03 are correct, an important consequence 
of the jet dominance in XTE J 11 18+480 is that all hard 
state sources are jet-dominated (in the sense that the jet 
power dominates over the X-ray power). This jet domi- 
nance also implies that all hard state sources should be 
radiatively inefficient. 

It is not clear whether most of the accretion energy is 
channeled into the jet or across the black hole event hori- 
zon. The L x oc m 2 j j dependence of the X-ray luminosity 



© 2002 RAS, MNRAS 000,0-?? 



Jet/disc coupling through energy reservoir 13 



(also observed in low- luminosity AGN, as inferred from 
the multivariate radio- X-rays-mass correlation, Merloni, 
Heinz & Di Matteo 2003; Falcke, Kording & Markoff 
2003) is similar to what is predicted by radiatively in- 
efficient accretion model. The reason for this inefficiency 
could be advection into the jet as well as advection into 
the black hole. 



large inner disc radius, implying a large total accretion 
rate (m ~ 0.2), which might be in conflict with the hy- 
pothesis of standard advection dominated flow theory. In 
the second case, in order to reproduce the very low in- 
ner disc temperature, an (uncomfortably) extreme coro- 
nal power is needed, together with substantial relativistic 
bulk motion of the coronal plasma, both possibly associ- 
ated with very high magnetic viscosity. 



6 CONCLUSIONS 

We have shown that the puzzling optical/X-ray correla- 
tions of XTE JIII8+480, can be understood in terms 
of a common energy reservoir for both the jet and the 
Comptonizing electrons. For illustration purpose, we have 
presented a specific shot noise variability model that re- 
produces all the main observed features of the multi- 
wavelength correlated variability. 

Our time dependent model is fairly general, and our 
main conclusions hold regardless of the specific geomet- 
rical and dynamical properties of the system. The main 
results can be summarized as follows: 

• Any energy reservoir model for XTE J1118+480 re- 
quires that the total jet power dominates over the X-ray 
luminosity. In particular, assuming that the compact jet 
synchrotron emission extends up to the optical band, this 
implies a radiative efficiency for the jet ej < 3 x 10~ 3 . Fol- 
lowing the same line of arguments as FGJ03, we showed 
that this situation is likely and probably represents a 
common feature of all black holes in the low-hard state 

• The range of typical variability timescales of dissipa- 
tion rate of the X-ray emitting plasma is consistent with 
the dynamical times of an accretion disc between ~10 
and 100 Schwarzschild radii (for a 10 Mq black hole). As 
expected, the X-ray variability can be associated with ei- 
ther a hot, thick inner accretion flow or with a patchy 
corona. 

• As jet launching requires large-scale coherent mag- 
netic structures, the energy dissipation in the jet should 
vary on longer time scales. Indeed, the derived range of 
timescales associated with the dissipation in the jet ex- 
tends to timescales more than 10 times larger than the 
X-ray emitting one. 

• By combining the information obtained from our 
time variability model with (time averaged) spectral anal- 
ysis we conclude that, whatever the accretion geometry, 
the whole disc-corona-jet system must be radiatively in- 
efficient. It is therefore possible, in principle, that the sys- 
tem is energetically dominated by the jet, as suggested by 
FGJ03, although whether the bulk of the accretion power 
is advected into the black hole or into the jet remains at 
present an open question. 

• In terms of specific dynamical models, we conclude 
that the observed properties of XTE J1118+480 during 
its low/hard state outburst are consistent with either 
an inner hot, quasi spherical, radiatively inefficient flow, 
from which the jet originates (Meier 2001), surrounded 
by a geometrically-thin, optically-thick cold disc, or with 
a powerful patchy, outflowing corona on top of an ex- 
tremely cold standard thin disc. In the first case, multi- 
colour disc fits to the UV/EUV spectrum indicate a very 
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